Three molecules of ubiquinone bind specifically to mitochondrial cytochrome bc1 complex.

Bifurcated electron flow to high potential "Rieske" iron-sulfur cluster and low potential heme b(L) is crucial for respiratory energy conservation by the cytochrome bc(1) complex. The chemistry of ubiquinol oxidation has to ensure the thermodynamically unfavorable electron transfer to heme b(L). To resolve a central controversy about the number of ubiquinol molecules involved in this reaction, we used high resolution magic-angle-spinning nuclear magnetic resonance experiments to show that two out of three n-decyl-ubiquinones bind at the ubiquinol oxidation center of the complex. This substantiates a proposed mechanism in which a charge transfer between a ubiquinol/ubiquinone pair explains the bifurcation of electron flow.

The pKa of His-24 in the folding transition state of apomyoglobin.

In native apomyoglobin, His-24 cannot be protonated, although at pH 4 the native protein forms a molten globule folding intermediate in which the histidine residues are readily protonated. The inability to protonate His-24 in the native protein dramatically affects the unfolding/refolding kinetics, as demonstrated by simulations for a simple model. Kinetic data for wild type and for a mutant lacking His-24 are analyzed. The pK(a) values of histidine residues in native apomyoglobin are known from earlier studies, and the average histidine pK(a) in the molten globule is determined from the pH dependence of the equilibrium between the native and molten globule forms. Analysis of the pH-dependent unfolding/refolding kinetics reveals that the average pK(a) of the histidine residues, including His-24, is closely similar in the folding transition state to the value found in the molten globule intermediate. Consequently, protonation of His-24 is not a barrier to refolding of the molten globule to the native protein. Instead, the normal pK(a) of His-24 in the transition state, coupled with its inaccessibility in the native state, promotes fast unfolding at low pH. The analysis of the wild-type results is confirmed and extended by using the wild-type parameters to fit the unfolding kinetics of a mutant lacking His-24.

Re-face stereospecificity of methylenetetrahydromethanopterin and methylenetetrahydrofolate dehydrogenases is predetermined by intrinsic properties of the substrate.

Four different dehydrogenases are known that catalyse the reversible dehydrogenation of N5,N10-methylenetetrahydromethanopterin (methylene-H4MPT) or N5,N10-methylenetetrahydrofolate (methylene-H4F) to the respective N5,N10-methenyl compounds. Sequence comparison indicates that the four enzymes are phylogenetically unrelated. They all catalyse the Re-face-stereospecific removal of the pro-R hydrogen atom of the coenzyme's methylene group. The Re-face stereospecificity is in contrast to the finding that in solution the pro-S hydrogen atom of methylene-H4MPT and of methylene-H4F is more reactive to heterolytic cleavage. For a better understanding we determined the conformations of methylene-H4MPT in solution and when enzyme-bound by using NMR spectroscopy and semiempirical quantum mechanical calculations. For the conformation free in solution we find an envelope conformation for the imidazolidine ring, with the flap at N10. The methylene pro-S C-H bond is anticlinal and the methylene pro-R C-H bond is synclinal to the lone electron pair of N10. Semiempirical quantum mechanical calculations of heats of formation of methylene-H4MPT and methylene-H4F indicate that changing this conformation into an activated one in which the pro-S C-H bond is antiperiplanar, resulting in the preformation of the leaving hydride, would require a deltadeltaH(f) of +53 kJ mol-1 for methylene-H4MPT and of +51 kJ mol-1 for methylene-H4F. This is almost twice the energy required to force the imidazolidine ring in the enzyme-bound conformation of methylene-H4MPT (+29 kJ mol-1) or of methylene-H4F (+35 kJ mol-1) into an activated conformation in which the pro-R hydrogen atom is antiperiplanar to the lone electron pair of N10. The much lower energy for pro-R hydrogen activation thus probably predetermines the Re-face stereospecificity of the four dehydrogenases. Results are also presented explaining why the chemical reduction of methenyl-H4MPT+ and methenyl-H4F+ with NaBD4 proceeds Si-face-specific, in contrast to the enzyme-catalysed reaction.

N-carboxymethanofuran (carbamate) formation from methanofuran and CO2 in methanogenic archaea. Thermodynamics and kinetics of the spontaneous reaction.

N-carboxymethanofuran (carbamate) formation from unprotonated methanofuran (MFR) and CO2 is the first reaction in the reduction of CO2 to methane in methanogenic archaea. The reaction proceeds spontaneously. We address here the question whether the rate of spontaneous carbamate formation is high enough to account for the observed rate of methanogenesis from CO2. The rates of carbamate formation (v1) and cleavage (v2) were determined under equilibrium conditions via 2D proton exchange NMR spectroscopy (EXSY). At pH 7.0 and 300 K the second order rate constant k1* of carbamate formation from 'MFR'(MFR + MFRH+) and 'CO2' (CO2 + H2CO3 + HCO3-+ CO32-) was found to be 7 M-1.s-1 (v1 = k1* ['MFR'] ['CO2']) while the pseudo first order rate constant k2* of carbamate cleavage was 12 s-1 (v2 = k2* [carbamate]). The equilibrium constant K* = k1*/k2* = [carbamate]/['MFR']['CO2'] was 0.6 M-1 at pH 7.0 corresponding to a free energy change DeltaG degrees ' of + 1.3 kJ.mol-1. The pH and temperature dependence of k1*, of k2* and of K* were determined. From the second order rate constant k1* it was calculated that under physiological conditions the rate of spontaneous carbamate formation is of the same order as the maximal rate of methane formation and as the rate of spontaneous CO2 formation from HCO3- in methanogenic archaea, the latter being important as CO2 is mainly present as HCO3- which has to be converted to CO2 before it can react with MFR. An enzyme catalyzed carbamate formation thus appears not to be required for methanogenesis from CO2. Consistent with this conclusion is our finding that the rate of carbamate formation was not enhanced by cell extracts of Methanosarcina barkeri and Methanobacterium thermoautotrophicum or by purified formylmethanofuran dehydrogenase which catalyzes the reduction of N-carboxymethanofuran to N-formylmethanofuran. From the concentrations of 'CO2' and of 'MFR' determined by 1D-NMR spectroscopy and the pKa of H2CO3 and of MFRH+ the concentrations of CO2 and of MFR were obtained, allowing to calculate k1 (v1 = k1 [MFR] [CO2]). The second order rate constant k1 was found to be approximately 1000 M-1 x s-1 at 300 K and pH values between 7.0 and 8. 0 which is in the order of k1 values determined for other carbamate forming reactions by stopped flow.

Protonation behavior of histidine 24 and histidine 119 in forming the pH 4 folding intermediate of apomyoglobin.

Heteronuclear NMR methods are used to study the protonation of histidine and aspartate residues in the acid-induced unfolding of recombinant sperm whale apomyoglobin. The results are combined with fluorescence and circular dichroism measurements of acid-induced unfolding of wild-type and double mutant (H24V/H119F) proteins. They are consistent with a simple model in which the failure to protonate a single buried histidine, H24, is largely responsible for the partial unfolding of native (N) wild-type apomyoglobin to the pH 4 folding intermediate (I). H24 is known to form an unusual interaction in which its side chain is buried and hydrogen-bonded to the side chain of H119. Two-dimensional 1H-15N heteronuclear NMR spectra indicate that H24 is present in the rare delta tautomeric form and remains neutral until N unfolds to I, while H119 becomes protonated before the N --> I reaction occurs. In the H24V/H119F double mutant, all histidines are protonated in N and the N --> I reaction occurs at lower pH. Therefore, the protonation of aspartate and/or glutamate residues must provide an additional driving force for the N to I reaction. Two-dimensional 1H-13C NMR experiments are used to measure the protonation of aspartates in selectively 13C-labeled apomyoglobin; the results indicate that none of the aspartate residues has a strongly depressed pKa in N, as would be expected if it forms a stabilizing salt bridge.

Interaction of minor groove ligands to an AAATT/AATTT site: correlation of thermodynamic characterization and solution structure.

A combination of circular dichroism spectroscopy, titration calorimetry, and optical melting has been used to investigate the association of the minor groove ligands netropsin and distamycin to the central A3T2 binding site of the DNA duplex d(CGCAAATTGGC).d(GCCAATTTGCG). For the complex with netropsin at 20 degrees C, a ligand/duplex stoichiometry of 1:1 was obtained with Kb approximately 4.3 x 10(7) M-1, delta Hb approximately -7.5 kcal mol-1, delta Sb approximately 9.3 cal K-1 mol-1, and delta Cp approximately 0. Previous NMR studies characterized the distamycin complex with A3T2 at saturation as a dimeric side-by-side complex. Consistent with this result, we found a ligand/duplex stoichiometry of 2:1. In the current study, the relative thermodynamic contributions of the two distamycin ligands in the formation of this side-by-side complex (2:1 Dst.A3T2) were evaluated and compared with the thermodynamic characteristics of netropsin binding. The association of the first distamycin molecule of the 2:1 Dst.A3T2 complex yielded the following thermodynamic profile: Kb approximately 3.1 x 10(7) M-1, delta Hb = -12.3 kcal mol-1, delta Sb = -8 cal K-1 mol-1, and delta Cp = -42 cal K-1 mol-1. The binding of the second distamycin molecule occurs with a lower Kb of approximately 3.3 x 10(6) M-1, a more favorable delta Hb of -18.8 kcal mol-1, a more unfavorable delta Sb of -34 cal K-1 mol-1, and a higher delta Cp of -196 cal K-1 mol-1. The latter term indicates an ordering of electrostricted and structural water molecules by the complexes. These results correlate well with the NMR titrations and are discussed in context of the solution structure of the 2:1 Dst.A3T2 complex.

Complexes of the minor groove of DNA.

An increasing number of high-resolution structures suggest that both the minor and major grooves of DNA can function as receptors for proteins and small molecules. In this review, we try to illustrate the diversity of small molecule ligands that are capable of specifically recognizing the minor groove of DNA. Complex formation results in varying degrees of conformational changes in both DNA and ligands. The discussion focuses on intermolecular interactions that contribute to binding affinity and specificity. There probably is no simple general recognition code that explains the binding specificity of minor-groove ligands. To understand DNA recognition by small molecules, characterization of the binding mode at near-atomic resolution must be combined with thermodynamic data on the energetics of ligand binding to short oligonucleotides.

Design of a G.C-specific DNA minor groove-binding peptide.

A four-ring tripeptide containing alternating imidazole and pyrrole carboxamides specifically binds six-base pair 5'-(A,T)GCGC(A,T)-3' sites in the minor groove of DNA. The designed peptide has a specificity completely reversed from that of the tripyrrole distamycin, which binds A,T sequences. Structural studies with nuclear magnetic resonance revealed that two peptides bound side-by-side and in an antiparallel orientation in the minor groove. Each of the four imidazoles in the 2:1 ligand-DNA complex recognized a specific guanine amino group in the GCGC core through a hydrogen bond. Targeting a designated four-base pair G.C tract by this synthetic ligand supports the generality of the 2:1 peptide-DNA motif for sequence-specific minor groove recognition of DNA.

Short peptide fragments derived from HMG-I/Y proteins bind specifically to the minor groove of DNA.

Short peptides derived from chromosomal proteins have previously been proposed to bind specifically to the minor groove of A,T-rich DNA [for a review, see M. E. A. Churchill and A. A. Travers (1991) Trends Biochem. Sci. 16, 92-97]. Using NMR spectroscopy, we investigated the DNA binding of SPRKSPRK, which is one such A,T-specific motif. Under the conditions studied SPRKSPRK interacts only nonspecifically with d(CGCAAAAAAGGC).d(GCCTTTTTTGCG). The peptides TPKRPRGRPKK, PRGRPKK, and PRGRP derived from the non-histone chromosomal protein HMG-I/Y, however, bind specifically to the central A,T sites of d(CGCAAATTTGCG)2 and d(CGCGAATTCGCG)2. 2D NOE measurements show that the RGR segment of each peptide is in contact with the minor groove. The arginine side chains and the peptide backbone are buried deep in the minor groove, in a fashion generally similar to the antibiotic netropsin. Under the same conditions the peptide PKGKP does not interact with the same oligonucleotide duplexes, indicating that the arginine guanidinium groups are major determinants of the A,T specificity.

Structural and dynamic characterization of the heterodimeric and homodimeric complexes of distamycin and 1-methylimidazole-2-carboxamide-netropsin bound to the minor groove of DNA.

NMR spectroscopy combined with molecular modeling was used to characterize a heterodimeric complex with Dst and 2-ImN bound in the minor groove of d(GCCTAACAAGG).d(CCTTGTTAGGC) (1:1:1 2-ImN.Dst.DNA complex). The imidazole-pyrrole-pyrrole ligand 2-ImN spans 5'-GTTA-3' of the TAACA.TGTTA binding site with the imidazole nitrogen specifically recognizing the guanine amino group. The Dst ligand lies along the 5'-AACA-3' sequence and complements the 2-ImN ligand in the formation of the antiparallel side-by-side heterodimeric complex. Titrations of the same site with Dst or 2-ImN alone yield homodimeric complexes (2:1 ligand.DNA) of lower stability than the 1:1:1 2-ImN.Dst.DNA complex. Dst and 2-ImN binding to d(CGCAAACTGGC).d(GCCAGTTTGCG) was also investigated. The 1:1:1 2-ImN.Dst.DNA complex is again the most stable complex with the AAACT.AGTTT site and is similar to the TAACA.TGTTA complex. No monomeric binding of either 2-ImN or Dst was observed to either site.

Antiparallel side-by-side dimeric motif for sequence-specific recognition in the minor groove of DNA by the designed peptide 1-methylimidazole-2-carboxamide netropsin.

The designed peptide 1-methylimidazole-2-carboxamide netropsin (2-ImN) binds specifically to the sequence 5'-TGACT-3'. Direct evidence from NMR spectroscopy is presented that this synthetic ligand binds DNA as a 2:1 complex, which reveals that the structure is an antiparallel dimer in the minor groove of DNA. This is in contrast to the 1:1 complexes usually seen with most crescent-shaped minor groove binding molecules targeted toward A+T-rich tracts but reminiscent of a dimeric motif found for distamycin at high concentrations. These results suggest that sequence-dependent groove width may play an important role in allowing an expanded set of DNA binding motifs for synthetic peptides.

Covalent modification of guanine bases in double-stranded DNA. The 1.2-A Z-DNA structure of d(CGCGCG) in the presence of CuCl2.

We have solved the single crystal structure to 1.2-A resolution of the Z-DNA sequence d(CGCGCG) soaked with copper(II) chloride. This structure allows us to elucidate the structural properties of copper in a model that mimics a physiologically relevant environment. A copper(II) cation was observed to form a covalent coordinate bond to N-7 of each guanine base along the hexamer duplex. The occurrence of copper bound at each site was dependent on the exposure of the bases and the packing of the hexamers in the crystal. The copper at the highest occupied site was observed to form a regular octahedral complex, with four water ligands in the equatorial plane and a fifth water along with N-7 of the purine base at the axial positions. All other copper complexes appear to be variations of this structure. By using the octahedral complex as the prototype for copper(II) binding to guanine bases in the Z-DNA crystal, model structures were built showing that duplex B-DNA can accommodate octahedral copper(II) complexes at the guanine bases as well as copper complexes bridged at adjacent guanine residues by a reactive dioxygen species. The increased susceptibility to oxidative DNA cleavage induced by copper(II) ions in solution of the bases located 5' to one or more adjacent guanine residues can thus be explained in terms of the cation and DNA structures described by these models.

Base-specific binding of copper(II) to Z-DNA. The 1.3-A single crystal structure of d(m5CGUAm5CG) in the presence of CuCl2.

The single crystal structure of d(m5CGUAm5CG) soaked with copper(II) chloride was solved to atomic (1.3 A) resolution to study the base specificity of copper binding to double-stranded DNA. In the present copper(II) chloride-soaked structure, four crystallographically unique copper(II) complexes were observed bound to five of the six purine bases in the hexamer duplex. Covalent copper(II) binding occurred at N-7 of all four guanine bases and at one of the two adenine bases in the DNA duplex. Copper binding was not observed at the position (Ade4) located in an open solvent channel, whereas the second adenine site (Ade10) shared a complex with a guanine residue (Gua12) of a neighboring symmetry-related hexamer. The coordination geometries and distribution of these copper(II) complexes at the guanine bases in the crystal were comparable to the analogous sites in the isomorphous copper(II) chloride-soaked d(CGCGCG) crystal (Kagawa, T., Geierstanger, B. H., Wang, A. H.-J., and Ho, P.S. (1991) J. Biol. Chem. 266, 20175-20184). Thus, the decreased copper(II) binding affinity for Ade4 was not an artifact of crystal packing, but is intrinsic to the chemical properties of this purine base in duplex DNA. This suggests that the adenine bases in dilute solutions of Z-DNA and more generally other duplex DNA conformations are not susceptible to copper(II) modification. Thus, preferential copper(II) binding at guanine bases over adenine bases in double-stranded DNA may explain the observed specificity of copper(II)-induced oxidative DNA damage near guanine residues (Yamamoto, K., and Kawanishi, S. (1989) J. Biol. Chem. 264, 15435-15440; Sagripanti, J.-L., and Kraemer, K. H. (1989) J. Biol. Chem. 264, 1729-1734). The sharing of a single copper(II) complex by Ade10 and Gua12 of an adjacent hexamer suggests that additional and perhaps specific DNA-DNA interactions, as may be found in the densely packed environment of the nuclear matrix in the cell, may render N-7 of adenine bases prone to copper(II) modification.